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Recent experiment on the tunneling spectra of ultrasmall superconducting grains revealed an 
unusual structure of the lowest differential conductance peak for grains in the odd charging states. 
We explain this behavior by nonequilibrium "gap less" excitations associated with different energy 
levels occupied by the unpaired electron. These excitations are generated by inelastic cotunneling. 
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The electron-electron interaction strongly affects the 
tunneling transport through ultrasmall metallic grains. 
Many phenomena of grains in the normal state are de- 
scribed by the orthodox model in which the interaction 
takes a simple form E^^^^ = Q'^/{2C), where Q and 
C are the total charge and capacitance of the system 
respectivelytl. This model is successful in describing the 
Coulomb blockade, which is essentially the quantization 
of the number of electrons in the grain, when the gate 
voltage is tuned away from the charge degeneracy point. 
Because of this quantization, the zero-bias conductance 
of the system vanishes, while the current / as the func- 
tion of the source-drain voltage V shows a threshold be- 
havior. Fine structure of the current -voltage cur|S{e was 
attributed to single-electron levels of the systemHa. 

However, there is a variety of effects which are not de- 
scribed by the orthodox model. They may become espe- 
cially pronounced when the system is driven out of equi- 
librium, thereby exploring many excited states. For nor- 
mal metallic grains, nonequilibrium steady states were 
shown to lead to clustering of resonances in the differ- 
ential conductance due to fluctuations of the interaction 
energy among electron^. The experimental findings for 
normal metallic grainsO can be summarized as follows: 

1. Clusters of resonances of the differential conduc- 
tance spectra appeared at high source-drain voltage 
for which nonequilibrium configurations of the elec- 
trons become energetically allowed. The first peak, 
however, did not split into several resonances. 

2. The width of each resonance cluster was found to 
be much smaller than the mean level spacing d. It 
is of order d/g, where 5 ~ 5 is the dimensionless 
conductance of the grain. 

Recently, Ralph et measured the tunneling reso- 
nance spectra of ultrasmall superconducting grains. The 
number of electrons in the system was controlled by a 
gate voltage. The results of this experiment show that: 

f . For the ground state of the grain with an even num- 
ber of electrons, the first peak of the differential 
conductance is merely shifted by the gate voltage 
Vg. The shape of this peak does not change over 



2. 



a large interval of Vg. Contrarily, if the grain con- 
tains an odd number of electrons, the height of the 
first peak rapidly reduces with a change of the gate 
voltage. Moreover, a structure of subresonances de- 
velops on the low-voltage shoulder of this peak, see 
Fig. g. 

The characteristic energy scale between subreso- 
nances of the first peak is of the order of mean 
level spacing d. 
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FIG. 1. A schematic illustration of the differential conduc- 
tance of an "odd" superconducting grain as function of the 
source-drain voltage V, at various gate voltages Vg. Higher 
resonances are separated by the superconducting gap from 
the first one, and subresonances are developed as the first 
resonance is shifted by the gate voltage. 

These observations contrast the results for the normal 
case, nevertheless, it was suggestecB that the substruc- 
ture of the first peak is still associated with nonequilib- 
rium steady states of the grain. The purpose of this paper 
is to clarify the origin of these nonequilibrium states and 
the mechanism which generates them. 

The principle difference between odd and even grains is 
that all excitations of the latter are of energy larger than 
the superconducting gap 2 A. Therefore, a source-drain 
voltage in the range V < A/e can not induce nonequilib- 
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rium states. Odd grains, on the other hand, contain one 
unpaired electron, which may be shifted to various sin- 
gle electron levels with characteristic energy scale smaller 
than the mean level spacing d. For this reason even small 
source -drain voltage d < eV < A is sufficient to ex- 
cite the grain. The mechanism of excitation is inelastic 
cotunnelingJ. Tunneling into the excited grain requires 
less energetic electrons, and lead in turn to the substruc- 
ture on the low-voltage shoulder of the of the first res- 
onance, see Fig. 1. A closely related problem was con- 
sidered by Averin and Nazarovcl, however, their theory 
assumed that relaxation processes prevent the formation 
of nonequilibrium states. For anall metallic grains re- 
laxation processes are very slo-wQ, and therefore will be 
neglected in our theory. 

To describe the effect quantitatively, we construct the 
master equations governing the time evolution of prob- 
abilities of different electronic configurations of super- 
conducting grains allowing for second order cotunneling 
processes. The solution of these equations for two lim- 
iting cases (one in which two levels participate in the 
transport, and the other when a large number of levels 
contribute) explains the substructure of the first peak of 
the differential conductance illustrated in Fig. 1. 

Our model Hamiltonian is given hy H = Hq + Ht + 
^int- Here Hq describes the noninteracting electrons in 
the left (L) and right (i?) leads and in the metallic grain. 
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Single particle energies ^^q, are measured relative to the 
chemical potentials of the left /il — aeV and right 
fiR = {a—l)eV leads (the numerical coefficient < a < 1 
depends on the capacitances between the grain and the 
leads). Tunneling across the barriers is described by 



where T^^J'' are the tunneling matrix elements. The in- 
teraction Hamiltonian is 



E E^ii^^<^^ + H-c. 

'=L,R q,a 
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Ant = ^(^-^o)^-^E4T^W 
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The first term in Eq. (g) is the orthodox model for the 
repulsive part of the electron-electron interaction: C is 
the total capacitance of the grain, N = CjCj is the 
number of electrons in the dot, and A^o is a continuous 
parameter controlled by the gate voltage. Nq determines 
the finite charging energy required to insert, {7+, or to 
remove, C/_, one electron. 



U±^-[^±{N-No)], |A^-Ao|<i. (4) 

The second term in Eq. (^ drives the grain to the super- 
conducting state (we assume zero temperature). 



We consider the experimentally relevant case, e^/C ^ 
A, so that the grain has well defined number of elec- 
trons. If this number is even N = 2m, the ground state 
energy (we will omit charging part of the energy and re- 
store it later), i?2m — ^2m +r3w"-i can be calculated in 
the mean field approximationElilj by minimizing thermo- 
dynamic potential — J2ki^k ~ ^k) + A^/A where 
ej. = (^2 + A2)i/2^ with respect to A {dF/dA = 0), and 
by fixing the chemical potential via 2m = —dF/dfx. All 
the excited states of even dots are separated from the 
ground state by a large energy, 2 A. Considering now the 
energy spectrum of an odd grain, N — 2m — 1, we notice 
that the second term in Eq. (^ operates only within spin 
singlet states. Therefore, to calculate the low-lying ex- 
cited states in this case, we fill the single-electron state j 
with one electron, and then find the ground state of the 
remaining 2to electrons with state j excluded from the 
Hilbert space. In the mean field approximation it cor- 
responds to the minimization of the thermodynamic po- 
tential i^2m-i = Efc5^j(6-efc) + ^V^ + Cj- The excited 
states with energies smaller than A are characterized by 

a single index, j and will be denoted by i?2m-i- -^'^ what 
follows we will need the energy cost of introducing an 
additional electron into the odd state: C/+ -I- Cj, where 
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^2m-i- Assuming A ^ d, one obtainsi 
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We turn now to the kinetics of a superconducting grain. 

2 

Consider the regime where U+ ~ U A, t/_ ^ ^ U, 

2 

and ^ ^ A 3> d. We also assume the conductance of 
the tunnel barriers to be much smaller than /h, and 
that the source-drain voltage is small eV < A. The sim- 
plicity brought to the problem in this regime of param- 
eters stems from the fact that there is only one avail- 
able state with an even number of electrons (because 
U- ^ U+ one can only add an electron to grain but not 
subtract one), and whenever the grain contains an even 
number of electrons it is in its ground state. This im- 
ply that even grains cannot be driven out of equilibrium 
state, while for odd grains tunneling (and cotunneling) 
takes place via unique state. 

Henceforth, we concentrate on grains with an odd 
ground state. Let us denote by the probability of 
finding the grain with an even number of electrons, and 
by Pj the probability to find the grain in the odd state 
j. Since these states are spin degenerate in the absence 
of magnetic field, Pj will denote the sum Pjj + Pj^i- 
The master equations for the probabilities PcPj have 
the form 
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dt 
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where ToX 



and VyXo are the transitions rates from the 
odd j-th state to the even and from the even to odd re- 
spectively, while Ti_,j is the rate of transition from the 
i-th to the j-th odd states. Equations (|^) are not lin- 
ear independent, so they have to be supplied with the 
normalization condition Pe + = 1- Current in the 

steady state equals to the electron flow through, say, the 
left barrier, and for positive V it is given by 
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Transition from the j-th odd state into the even state 
occurs when fj.^ > U + Ej. The amplitude of this tran- 
sition is calculated by first order perturbation theory in 
the tunneling Hamiltonian (0). Fermi's golden rule yields 



'^e = 9L 
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where (Jl is the dimensionless conductance of the left 
tunnel barrier per one spin, Uj = (1 + /ej)/2 is the 
coherence factor, 9{x) is the unit step function, and 
PLj = fl\ipjirL)\'^ , where fl is the volume of the grain 
and i^j{rL) is the value of j-th single particle wave func- 
tion at the left point contact ■ Energies are given by 
Eq. (Q) and U = U+ is defined in Eq. (||). Similarly, the 
rate of transition from even state to i-th odd state, by 
tunneling of an electron from the dot to the right lead, 
is given by 
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where is the dimensionless conductance of the right 
tunnel barrier, Vi = {I - ^i/e.i)/2, and = ri|V'i(r-i?,)P, 
where rn is the position of the right point contact. 

A change in the occupation confisuration of the odd 
states occurs via inelastic cotunnclingZI. This mechanism 
is a virtual process in which an electron tunnels into j-th 
available level and another electron tunnels out from the 
i-th level. Calculating this rate by second order pertur- 
bation theory in the tunneling Hamiltonian, one obtains 



for eV > e.i 



gLgRCpU^vf pLjPBi(eV - £j + £i) 

Stt^MU + £j - pl){U + e, - ^^R) 
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Ei, pl < U + Ej, pn < U + Ei, and zero 
otherwise. Tj^i diverges in the limits fiL ^ U + Ej and 
pn ^ U + Ei. It signals that a real transition takes over 
the virtual one. The region of applicability of Eq. (p^) 
is, therefore, U + Ej — p,L > 7 and U + Ei — > 7 where 
7 ^ gd/An is the width of a single particle level in the dot 
due to the coupling to the leads, g = gL + 9r- However, 
the interval of biases where Eq. ( p^ ) is not valid is nar- 
row, and to the leading approximation in Ht our results 
will be independent of this broadening. 
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FIG. 2. Inelastic cotunneling process can drive an "odd" 
superconducting grain out of its ground state. In the ground 
state, the single particle level indicated by [/ -I- £0 is occu- 
pied by one electron. Excited states are those in which the 
unpaired electron is shifted to other single particle levels. In 
a nonequilibrium steady state, low single particle levels be- 
come available for resonant tunneling, leading to a subres- 
onances structure of the differential conductance shown in 
Fig. 1. State j shown to be filled with two electrons should 
be understood as a coherent superposition of double occupied 
and empty states with weights and v?j respectively. 

Let us now apply Eqs. (^) and (0) to describe the ap- 
pearance of the low- voltage substructure of the first peak. 
We will consider two situations: (i) small voltage such 
that only one subresonance can emerge on the shoulder 
of the leading one, and (ii) large voltage, d <C eV < A, 
where the substructure of the main resonance consists of 
a large number of subresonances. 

In the first case, the chemical potentials of the left and 
right leads are such that transport through the grain in- 
volves only two levels: £0 and Ei < Eq corresponding to 
the ground and the first excited states of the odd grain. 
We solve Eqs. (^) for probabilities Po,Pi and Pe using 
Eqs. (|^) and (^). There are two distinct regimes of the 
source-drain voltage: (1) fiL < U + Eo where transport is 
dominated by cotunneling, and (2) pl 1^ U + Eq where 
state "0" is available for resonant tunneling. The sub- 
structure of the first resonance in the differential conduc- 
tance appears in the first regime. Below we show that as 
PL passes through U + Ei, see Fig. ||, there is a discon- 
tinuity in the current-voltage curve. In the first regime, 
the total current to the leading approximation in gh^dR 
is a sum of two contributions, / ~ leq -f Ine- The first, 
leq — cFo-jo + 2ero^i, is the equilibrium current coming 
from cotunneling. The second contribution is associated 
with the nonequilibrium population of state "I" and is 
given by 
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Assuming that the voltage drop eV ~ pL ~ Pr is larger 
than the energy difference d = eq — Ei, the jump in the 
nonequilibrium current is: 
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sine = Cie- 
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where ci = 4,UqViPloPri is a constant of order unity. 
This jump in the nonequihbrium current leads to the 
peak in the differential conductance spectra. Formula (|p) 
has simple interpretation. Up to numerical prefactors it 
is a product of two factors: first is the probability of find- 
ing the grain with an unpaired electron in state "1". ft 
is proportional to gii{d/d){l — d/eV), and increases with 
the voltage V and as d = Eq — £i ^ 0. The second factor 
is associated with the rate in which the state "1" is filled 
with an electron, egLd/h. 

The magnitude of the jump (||) should be compared to 
the jump in the current as /i^ increases above U + eo, and 
real transition via the even state become allowed. To the 
leading order in and gn, the current in this regime is 
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where C2 is a constant of order unity having structure 
similar to ci. Comparing the current jump, Sine, with 
that associated with the resonant tunneling, 51, we find 
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Thus nonequihbrium population of the excited level of 
the odd-grain leads to the appearance of a subresonance 
at small V, however, its height is much smaller than that 
of the main resonance. 

We turn now to the second regime of the parameters, 
d <C eV < A, in which many levels contribute to the 
transport. Again, we focus our attention on the cotun- 
neling regime, p^ < U + eo. We show that the character- 
istic amplitude of the subresonances in this regime may 
become comparable to the amplitude of the main peak. 

To the leading order in 5l,5_r, and d/A, the steady 
state solution of the rate equations at /i^ = [/ -I- ei -|- is 
Po — 1, while for the other probabilities we have 
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The characteristic number of states contributing to the 
current (^ is large as y/AeV /d so that mesoscopic fluc- 
tuations of the tunneling rates and of the inter- level spac- 
ings may be neglected. Additional large factor, V AeV /d, 
comes from the summation over the levels in Eq. (12), 
and we find 



PL = U + Ej+Q 



(13) 



Once again, the current jumps each time pL passes 
through U + ej. This jump for large j (but still such 
that U + Ej — PL ^ eV) scales as and the ratio of 



the jump at j = 1 to the jump at the resonance level ( [lO| ) 
is given by 
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, d<t:eV<A. (14) 



Noticing that d = Eq — £i ~ d"^ /2A, we see that the first 
subresonance becomes comparable in height to the main 
one at voltages as small as eV ~ 4TT^d^/A'^{gL + iga). 

We conclude by comparing the above results with the 
experimental data of Ref. [ g] . There A « 5d, the con- 
ductances in the normal state are gn ~ 10(?l ~ 1/8, and 
the leads are also superconducting. The singularity in 
the density of states of the leads imply that the effective 
conductance is increased by factor of 2 — 3. Neglecting 
inelastic processes and the Josephson coupling, these pa- 
rameters imply that when eV « 2d the ratio of the sub- 
resonances amplitude to that of the first resonance is of 
order one, while at eV ^2d~ d/5 it is of order of 1%. It 
implies that first subresonance peak associated with tun- 
neling into state "1" cannot be resolved, both because its 
amplitude and its distance from the main peak are too 
small. However next subresonance appear already at dis- 
tance of order d from the main resonance, and for V > 2d 
have an amplitude comparable with the main resonance. 
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